Quantum limit of photothermal cooling 
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We study the problem of cooling a mechanical oscillator using the photothermal (bolometric) 
force. Contrary to previous attempts to model this system, we take into account the noise effects 
due to the granular nature of photon absorption. This allows us to tackle the cooling problem down 
to the noise dominated regime and to find reasonable estimates for the lowest achievable phonon 
occupation in the cantilever. 



I. INTRODUCTION 



The optomechanics of deformable cavities was born in 
the seventies from the seminal works of Braginsky 
In recent years the first successful observation of self- 
cooling due to photothermal @ and radiation pressure 
forces started a race to reach the quantum regime 
of mechanical motion 

The self-cooling process usually involves an optical cav- 
ity whose mirrors feel an optically-induced force propor- 
tional to the cavity population. This force depends on the 
mirror position because the cavity population depends on 
the cavity length. If this dependence is delayed in time, 
e.g., because the force depends on the position of the 
mirror at an earlier time, then self-cooling of the mir- 
ror motion can be achieved (given the right parameter 
regimes [HI]). 

These optomechanical self-cooling schemes can be clas- 
sified into different categories according to the nature of 
the optically-induced force that dominates the cooling 
process. The most commonly used is the radiation pres- 
sure force that exploits the pressure exerted by photons 
bouncing off the mirror. In this case, the delay in the 
force is given by the cavity storage time; that is, the 
time taken by the photon population inside the cavity to 
adjust when the cavity length is modified. The second 
kind of optically-induced force exploits the thermal de- 
formation of the mirror due to the absorption of photons 
19]. In this case, the delay is given by the heat dif- 
fusion time through the mirror. This force, dubbed the 
photothermal (or bolometric) force, has been observed in 
multiple experiments and can be much stronger than the 
radiation-pressure force in appositely designed systems 

0.HJE3. 

When either cooling mechanism brings the system to- 
wards the quantum regime, it also becomes necessary to 
consider (together with the dynamically-induced cooling 
rate) the spontaneous heating rate given by the intrinsic 
quantum noise. This is because photons are reflected, 
or absorbed, one-by-one in both cooling schemes. In 
the case of radiation-pressure-cooling, a quantum theory 
exists [22h26| ] showing that, with the ideal parameters, 
cooling to the ground state of the mirror motion is still 
possible. However, unlike the radiation pressure force, 
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FIG. 1: (Color online) Schematic diagram of an optomechan- 
ical system. The mechanical system, or cantilever, is fixed 
on the (blue) left edge, and operates as a movable mirror 
(right side, in orange) forming one half of an optical cavity 
(shown by the red constriction). The grey disk is a non- 
moveable (i.e., fixed) half-mirror. The optical and mechanical 
modes are coupled via both radiation-pressure and photother- 
mal forces. The optical cavity (bounded by the fixed grey disk 
and the (orange) right edge of the cantilever) has equilibrium 
frequency cj c , quality factor Q c , equilibrium length L c , photon 
lifetime F c , and is pumped by an external laser with power 
P and frequency ui p . The cantilever has thermal conductiv- 
ity «, length L m , and surface area s. The cantilever mode 
we consider for cooling has frequency u m , effective mass m, 
and quality factor Q m . The right edge of the cantilever ab- 
sorbs photons at the rate a, and has a thermal deformation 
coefficient \- 



until now the photothermal cooling process has not been 
investigated in the quantum limit. This is due to the in- 
trinsic difficulty in building a fully quantum theory of a 
process that is fundamentally dissipative and involves a 
macroscopic number of degrees of freedom (e.g., photon 
absorption, heat diffusion, and thermal deformation). 

Our goal in this paper is to construct a theory (albeit 
a phenomenological one) rigorous enough to be able to 
firmly answer questions of practical importance for ex- 
periments that employ photothermal cooling. Mainly we 
want to know whether it is possible to exploit photother- 
mal cooling to reach the quantum regime, and to identify 
the bounds on parameters necessary to reach this regime. 
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The theory we construct uses the model of a cantilever- 
mounted mirror [see Fig.[T] in order to be able to compare 
with specific experiments. However, in principle one can 
also apply our arguments to other systems, like microdisk 
optomechanical systems [4(J 5l[ . 

This paper is structured as follows: in Section [TT] we 
lay down our theory and analyze the results. After hav- 
ing described the issue of residual heating due to photon 
absorption in Section Um we will test our theory with nu- 
merical data taken from a well documented experiment 
in Section HVl 



II. THEORY OF PHOTOTHERMAL COOLING 
A. Basic approach 

Our aim is to describe the process of self cooling of an 
oscillation mode of a cantilever (see Fig.[IJ. In the follow- 
ing, unless explicitly stated, we will always be considering 
this specific mode. It is usually the fundamental (lowest 
frequency) mode of the cantilever, but this depends on 
the geometry of the optical cavity and the cantilever [2l[ . 

In the specific model we consider here, the cantilever 
acts as one of the mirrors of a pumped photonic cavity. 
The cantilever mode couples to the photonic field of the 
cavity because the frequency of the cavity, and thus its 
equilibrium population, depend on the cantilever posi- 
tion. Assuming the pumping laser to be red detuned, 
and referring to the coordinate axes as shown in Fig. [TJ 
we have that a shift of the cantilever toward positive x 
will cause a lowering of the cavity frequency and thus an 
increase in the photonic population (and vice versa). 

Even with a very reflective coating there will always 
be a residual photon absorption in the cantilever. When 
a photon is absorbed in the cantilever it will create an 
excess of heat (and thus an excess phonon population) 
in the region of the cantilever that serves as the optical 
cavity mirror. This heat will then diffuse through the 
cantilever following normal rules of heat transfer. 

Anharmonicity in the interatomic potential in the can- 
tilever give s rise to the phenomenon of thermal defor- 
mation [271 ] . As the average energies in the vibrational 
degrees of freedom of the atoms (and thus the phonon 
population) increase, the average interatomic distances 
change. If this deformation is not uniform through the 
whole cantilever, (e.g. if the cantilever is formed by layers 
of materials with different thermal expansion coefficients) 
then strains develop that exert an effective force on the 
entire cantilever. 

The number of absorbed photons, and thus the cor- 
responding excess heat, is proportional to the total 
photon population of the optical cavity. However the 
deformation-generated-force is exerted only after the 
heat has diffused through the cantilever. This implies 
that the cantilever is subject to a force dependent on its 
position at past times. 

The cavity photons act on the cantilever also through 



the well known radiation pressure effect, that is, they ex- 
ert a force on the cantilever as they bounce on the mov- 
able mirror. In usual radiation-pressure cooling experi- 
ments the cavity response time is of the order of the can- 
tilever mode frequency (the good cavity regime [231 l2~4l] ) . 
that is, the time the cavity population takes to adjust 
to the new cantilever position is of the same order of the 
cantilever period. This also results in a delayed force act- 
ing on the cantilever. However, in the systems optimized 
for photothermal cooling which we treat in this paper, 
the cavity response time is much shorter than the can- 
tilever oscillation period (the bad cavity limit) and we 
can thus safely neglect the retarded nature of such force. 

As mentioned in the introduction, to develop a mi- 
croscopic quantum theory of the whole process is a 
formidable task, due to the intrinsic many-body and dis- 
sipative nature of the phenomena involved (photon ab- 
sorption and heat transfer). However, given the time 
scales of the cantilever motion and of the thermal de- 
lay, we can assume that the quantized cantilever mode 
is evolving under the action of a semi-classical delayed 
feedback force. In other words, the whole system (with 
thermal delay) effectively measures the cantilever posi- 
tion (via the number of absorbed photons) , stores it in a 
classical signal (the heat) and transmits it to a classical 
actuator (the cantilever deformation). 

The cooling effect thus arises through a semi-classical 
feedback or backaction effect, not different in principle 
from electronically controlled active classical feedback 
mechanisms used in quantum optics [191, |28H30j . In order 
to better understand this parallelism it is worthwhile to 
look at the orders of magnitude involved in usual pho- 
tothermal experiments. For the cantilever in Ref. [20 the 
cantilever is 220 /im long and its thermal diffusion time 
is of the order of 0.5 ms. Comparing with the lifetime 
of phonons in silicon, which are in the picosecond scale 
[3 If , it is easy to understand that no quantum coherence 
survives through the effective feedback loop. 

Here, we thus construct a Langevin description of the 
cantilever under the action of the radiation pressure force 
and of the time-delayed (effective feedback) photother- 
mal force with an added noise source in the effective feed- 
back loop due to the photon absorption shot noise. 

B. Treatment of the photothermal force 

It is well known that the photothermal force, in the 
classical regime, can be modelled using a time-delayed 
force, depending on the position of the cantilever at past 
times. This approach has been proven to describe well 
experiments not only in the cooling regime 0, [20j | but 
also in the regime where self-sustained oscillations de- 
velop [34| . The force can thus be written as a time inte- 
gral over a memory kernel of the instantaneous force due 
to the thermal deformation 

rp mm f, s eLFdef (t,x(t)) 

F P h[x\(t) = j^dr h(t - t) — , (1) 
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where the memory kernel associated to the thermal dif- 
fusion delay r is taken to be 



h(t) = [l-exp(-t/T)}e(t), 



(2) 



with Q(t) the Heaviside function. 

From Eqs. ([T]) and ^ we see that the instantaneous 
force F dc { coincides with the actual photothermal force 
F p h we would have for a vanishing delay time r 



F dei (t,x) 



lim F ph [x](t). 

T— »0 



(3) 



As shown in Refs. |32j and l33l. in a coherently pumped 
optical system with a closed feedback loop, it is suffi- 
cient to use a semiclassical theory in order to describe 
the photon absorption shot noise. We limit ourselves to 
a semi-classical description of the electromagnetic field 
in a coherent state. The easiest way to do so is to start 
from a quantum Hamiltonian for the pumped photonic 
cavity [2J] 

H c = huj c (\ - a) a + Ee iu ^a + Ee^'W, (4) 

where a is the annihilation operator for a cavity photon, 
lo c is the cavity frequency, L c its equilibrium length, uj p is 
the pump frequency, E is a pump term and we assume the 
loss rate of cavity photon to be given by T c . Substituting 
c-numbers for the photonic operators (d — > a), as usual 
for coherent states, allows us to calculate the population 
of the cavity. Given the large discrepancy between the 
frequency of the cantilever uj rn (usually in the kHz to 
MHz range) and of the photonic cavity response time r c 
(usually in the GHz to THz range) we are not going to 
consider explicitly the dynamics of the photonic mode, 
but we will assume that it responds instantaneously to 
any change in the position of the cantilever. 

We thus abtain the instantaneous photon population 
of the cavity as a function of the cantilever position x, 



n c (x) 



E 2 



[Wc(l 



L c ■ 



•ri/4- 



(5) 



Supposing that we have a well defined cavity frequency 
lo c , we can find the relation between the pump term E 
and the input power P by equating the number of in- 
coming photons P/hujp with the number of photons lost 
n c T c . We thus obtain 



E = 



IT C P 
Ahuj„ 



(6) 



where r c is the cavity photon lifetime. 

In order to proceed in the derivation of a quantum 
Langevin equation for the cantilever, which takes into 
consideration photon shot noise, we need to fix the de- 
pendency of the photothermal force upon the intensity of 
the cavity field n c (x). 

While the microscopic derivation of such a force would 
be extremely complex and sample-dependent, we expect 



it to be a function of the heat absorbed by the can- 
tilever. In other words, the instantaneous photon energy 
cj c (1 — times the current of absorbed photons I c (x, t). 
Here we will limit ourselves to the simple case of a lin- 
ear dependence that, as we will see, gives good results 
when compared with experiments. We will thus write 
the instantaneous deformation force as, 



F de( (t,x) = \frhJc 



1- — )I c (x,t), 



(7) 



where x is a phenomenological deformation coefficient of 
the cantilever with dimension of the inverse of a velocity. 

We note that the current of absorbed photons I c can be 
written as an average component and a fluctuating com- 
ponent due to shot noise, the average being proportional 
to the number of photons present in the cavity (itself 
a function of x) and the noise to a coefficient times a 
stochastic noise term r](t) [H, H3 



Ic(x,t) 



an c (x) 
an c (x) 



SI c (t,x) 
N(x) V (t), 



(8) 



where a is the absorption rate of the cantilever. 

In order to have a cooling effect the noise term in 
Eq. (|8]) needs to be small compared with the average 
value. Assuming the cantilever oscillations to be much 
smaller than the equilibrium cavity length L Cl we can 
thus insert Eq. ([5]) into Eq. ([7]) and make a first order 
expansion both x and ?y. We thus obtain 



F dei (t,x) ~ F + (VF)x + N r)(t). 



(9) 



As the different absorption events are completely un- 
corrected we can, in a time interval short enough for 
neglecting the variation of x, consider the number of ab- 
sorptions as a Poisson process. We thus define 



n p (t)= / dt'I c {x,t'), 
Jo 



(10) 



as the number of photons absorbed up to time t. Its 
momenta are given by 



{n p {tf) 



tan c (x), 
t 2 a 2 n c (x) 2 



(11) 

dt' / dt"{ V (t')ri(t")). 



Imposing that n p has Poissonian, time independent 
statistics over the time interval we consider here, we have 



(n p (t)) = (n p (t) 2 ) - (np(t)) s 



(12) 



We thus arrive at the following condition for the correla- 
tor of current fluctuation in the linear approximation 



A^i'MO) = an c (0)S(t'-t" 



(13) 



The current I c is thus characterized by a white noise 
spectrum (as in the semi-classical treatment of optical 
feedback [32J, |33j ) and can be written as 

I c (x,t) — an c (x) + y/an c (0) V (t), (14) 
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where the white noise term rj(t) has the correlator 

( V (t') V (t")) = S(t'-t"). (15) 

Putting together Eqs. (IT4l .([5|) and (O and applying 
basic analysis we find the following values for the coeffi- 
cients VF and N 



VF 



2 (0)axhuj c 2lo c A - A 2 - if 



A 2 



r? 



(16) 



where A = uj c — uj p is the cavity detuning. In section IID 
we use these terms directly in the equation of motion for 
the cantilever. 



C. Treatment of the radiation pressure force 

Radiation pressure force lends itself to a microscopic 
quantum treatment and it has been studied in various 
publications [23, 24] . These works show that the effect of 
radiation pressure cooling can be easily modelled by an 
effective damping rate r rp and an associated noise term 

The amplitude of r rp , that is the cooling capacity of 
the radiation pressure force depends on the ratio between 
the cavity damping and the cantilever frequency and thus 
becomes negligible in the bad cavity regime we consider 
here. The same is not true for the noise term, that in- 
stead only depends on the intrinsic magnitude of the force 
and thus leads to a nonnegligible equilibrium noise pop- 
ulation. 

From Ref. Hi] we find the following formulas for the 
optical damping term and for the equilibrium noise pop- 
ulation nU 



r rp 



4n c (0)fir c w 2 



A 



A 2 



iLl 
El 

4 



(A 2 



El)?' 

4 / 



4w„,A 



(17) 



Having the damping rate F rp and the noise population 
in this form, and because we have assumed the bad- 
cavity limit, means that we do not have to explicitly con- 
sider the dynamics of the optical cavity. 



sure force 
P 

x = — , 
m 

p = -mu? m x - (r m + T rp )p + F th + F ph [a;] 



(18) 
Fpj 



where T m and F t h are the dissipation and the fluctuation 
terms given by the coupling with the bath [35T - I371 ] . 

Using the expression for the linearized time-delayed 
force from Eq. §§§ in Eq. (fT8"f we obtain the following 
second-order dynamical equation 



x(t) + u? m x(t) + (r m - 

dt'h(t - 1') 



r Ip )x(t) = 

VF d . , 

171 dt 



m 
N d 
m dt 



(19) 



This is an integro-diffcrential equation, and the integral 
has a convolution. It thus becomes an algebraic equation 
in Fourier space 



2 2 



F th (u) F ip {oj) 



m(l 

Nr)(u)) 



x{uj) 



m m(l + iojt) 



(20) 



As we see from Eq. (f2"01 , the photothermal force results 
in a frequency-dependent damping and a renormalization 
of the cantilever frequency (see also the discussion in Ref. 
20l ). In the usual (rigid cantilever) regime in which the 
frequency of the cantilever, renormalized by the constant 
{uj = 0) part of the photothermal force, 



VF 

777 ' 



(21) 



does not differ much from w m , the essential part of the 
dynamics takes place at w = u m , and we can thus take 
the photothermal-induced damping to be 



r p h — 



tVF 



m(i 



J' 



(22) 



Notice that, as explained in Ref. l2fj, the value at which 
Eq. (|21[) vanishes corresponds to the onset of mirror in- 
stability and sets an upper limit to the strength of the 
photothermal force 



VF 



< 1. 



(23) 



D. System dynamics: feedback in a quantum 
Langevin equation 

Using the results from the previous sections, we can 
describe the dynamics of the system with a set of cou- 
pled Langevin equations describing the cantilever mode 
coupled to a thermal bath and evolving under the con- 
junct effect of the photothermal and the radiation pres- 



Neglecting quantum (zero-point) fluctuations in the 
photothermal force we can thus calculate the phonon 
equilibrium population due to the photothermal force 
alone (i.e., neglecting both the coupling to the thermal 
bath and the radiation pressure force in Eq. ([2U]) ) as 



,ph 



N 2 



2huj m TVF 



(24) 



In principle we can also derive the total equilibrium 
phononic population due to the joint effect of the ther- 
mal bath, the radiation pressure and the photothermal 
effect by integrating Eq. (|20p (under the rigid-cantilever 
approximation). However it is easier to calculate this 
total population via a simple rate equation, which gives 



V r? th 

tot _ L mi'm 



ph 



r, 



r rp + r 



P h 



(25) 



Reference 




w m 


Qm 


n th 

""m 


n^ m,n , Eq. ([30J 


Ref. 38 
Ref. 20 
Ref. 39 


275 /im 
220 urn 
3.9 fim 


6.5 MHz 
46 KHz 
3.4 MHz 


1.5 x 10 b 
2.2 x 10 3 
2.9 x 10 3 


1.2 x 10 b 
1.7 x 10 8 
2.2 x 10 6 


5 

6.4 x 10 5 
7.8 x 10 3 



TABLE I: Length L m , frequency ui m , quality factor Q m , ther- 
mal equilibrium population at 77 K and final population 
n m' mm from Eq. (|30[) for different cantilevers reported in the 
literature. 



wher 



,th ; 



is the equilibrium population at temperature 



T that, given the low frequencies of the cantilever modes, 
can be approximated as kT /hw m . 

Inserting Eqs. CEf)]), dTTJ) and (JUJ), into Eq. ^ and 
remembering that we are working in a regime where 
<C T p h, we can rewrite Eq. (|25t as 



r rp « r 



,N 



where 



kTmujl 1 + o4t 2 
hu) m Q m VF w m T 



(26) 



(27) 



is the classical population due to photothermal cooling 
and 



,N 



1 I+uLt 1 



a x^cL c u m r 2u r A - A 2 - S 



1 



2w m r 2u> r A - A 2 



4 



(28) 



is the population due to noise effects. The first line in 
Eq. (128[) stems from the radiation pressure noise con- 
tribution, and the second from the photothermal noise 
contribution. 



E. Analysis of the results: the quantum regime 



From Eq. (|27|) we can, using Eq. (123)) and optimizing 
over V-F, find (consistent with previously known results 
[U H3]) that the efficiency of classical photothermal cool- 
ing is maximal for Lo m r = 1 and its ultimate value de- 
pends on the quality factor of the cantilever 



r m. 



(29) 



The optimization of Eq. (|27[) gives us a lower bound on 
the phonon population 



n£' min > 



fy^rnQm 



(30) 



Thus our estimate for the minimal classical population 
for the cantilever is proportional to its thermal, equilib- 
rium occupation kT/huj rn divided by the quality factor 
of the mechanical oscillator Q m . Given the rather low 



frequencies of the mechanical modes involved, the condi- 
tion of having an average occupation number lower than 
one puts a rather harsh requirement on the mechanical 
quality factor. In Table I we compare some data taken 
from experimentally realized cantilevers. We see that 
for the stressed silicon cantilever studied in Ref. [38| the 
theory predicts, assuming optimal parameters and liquid 
Nitrogen temperatures, a thermal population almost in 
the quantum regime. Note that only the cantilever in 
Ref. [20 has actually been used for photothermal cool- 
ing, and thus it is the only cantilever for which we can 
estimate the deformation coefficient \ an d the noise con- 
tribution to the final population (see Section HV)) . 

Assuming that we have a cantilever quality factor Q m 
large enough to bring n^ n into the quantum regime, 
we are thus confronted with the noise contributions in 
Eq. (|2"5)l . The noise in Eq. (|2"5| depends on various pa- 
rameters, and the two noise contributions (the radiation 
pressure in the first line and the photothermal in the sec- 
ond) often have inverse dependencies upon the parame- 
ters. Thus naively minimizing one noise term can enlarge 
the other. We are thus obliged to carefully optimize the 
parameters in order to calculate the minimal noise pop- 
ulation. We start by rewriting Eq. (|28p as 



r c 1 + ujI 



-QlA- 



1 

2A 



1. 



(A + 7) 



(31) 



2Q C A -A z - - 



where we have defined 
.4 



LU m T 



X^cLc 

and introduced the renormalized detuning 

A = A/r c , 

and the cavity quality factor 

Qc = Wc/r c - 



(32) 



(33) 



(34) 



Again the first term in the numerator of Eq. pip is from 
the radiation pressure noise, and the second is from the 
photothermal noise. Thus, the new parameter A encap- 
sulates the way in which the two noise terms are inversely 
proportional to one another. Optimizing Eq. (1311) over A 
we obtain 




(35) 



6 



where the optimal value of A is given by 



.4 



opt 




a 



2T c Ql 1 



(36) 



This optimal value of A makes the two noise terms equal, 
and minimum. The second factor in Eq. (1351) can be 
further optimized over Q c , yielding its minimum for large 



T c l + u;2 T 2 



1 

4A' 



(37) 



Choosing a large enough detuning A we can thus ignore 
the second factor in Eq. (|37[) and we obtain our final 
estimate for the minimum noise population 



> 



a 2o;2 r 2 



(38) 



We see that, given a large enough cavity quality factor 
Qc and the optimal value for x from Eq. ([36]), two fac- 
tors determine the minimum noise population: the delay 
parameter w m r and T c /a. The second one, T c /a, is the 
ratio between the total photon loss rate (r c , the sum 
of the mirror absorbtion a and other radiative losses), 
and the loss part due to absorption in the moving mirror 
alone (thus this ratio is always larger than one). Equa- 
tion (|38|) tells us that, due to the interplay between the 
radiation pressure and photothermal noises, while it is 
possible to reach the "quantum regime" with appropri- 
ate parameters, it is not possible to get arbitrarily close 
to the ground state as in the case of radiation-pressure- 
dominated experiments |24j . In order to reach the quan- 
tum regime it is necessary not only to limit all the losses 
other than the cantilever absorption in order to lower 
r c /a, but also to design samples with rather large val- 
ues of Lo m T. Furthermore, while large values of w m T help 
lower the noise population, Eq. (|27|) tells us that they cor- 
respondingly reduce the classical cooling efficiency, thus 
increasing the classical phonon population. 

To reiterate, to reach the quantum regime we require 
the following conditions: 



• The feedback cooling term, Eq. (|30j) . tells us that 
we need a large cantilever quality factor Q m , large 
cantilever frequency w m , and low initial tempera- 
ture T, so that Q m » kT/huj m . 

• To simultaneously minimize the photothermal and 
radiation pressure noises, we need a high quality 
optical cavity Q c 3> 1, a large enough detuning A, 
and a deformation constant \ given by Eq. 



• These conditions give us Eq. (|38j) . the final equilib- 
rium phonon population due to noise, which is then 
minimized by choosing T c /a small. Again, this can 
further be minimized by increasing w m r, at the cost 
of reducing the feedback cooling efficiency, and thus 
requiring a higher cantilever quality factor Q m . 



III. RESIDUAL HEATING FROM PHOTON 
ABSORPTION 

Since the photothermal force is due to photon absorp- 
tion, the same process that gives the cooling is simulta- 
neously heating up the cantilever as a whole. This is a 
rather intriguing aspect of photothermal cooling that we 
have, until now, neglected. 

This heating mechanism does not modify the funda- 
mental bounds we derived, accounting simply for an in- 
creased temperature in Eq. (|27|) However, the photother- 
mal heating effect is important if we want to quantita- 
tively fit our model with experimental data, as it can 
noticeably influence the observed temperature. 

The exact description of such residual heating on the 
cantilever motion would depend on the microscopic struc- 
ture of the cantilever and it is thus not in the scope 
of this paper. Here we will limit ourselves to a rather 
phenomenological model, able to mimic the interesting 
physics, and depending on only one adjustable param- 
eter. The simplest way to take into account this phe- 
nomenon is to use, inside Eq. (|27j) . not the environmental 
temperature T, but an average equilibrium temperature 
under constant illumination (average because the tem- 
perature profile will in general vary over the mechanical 
mode region). 

We now consider a rectangular cantilever, character- 
ized by a constant thermal conductivity (constant in 
space and in temperature), ignoring radiative processes 
and assuming that the absorption processes happen ho- 
mogeneously over one of the surfaces (see Fig. [T] for a 
schematic illustration). Calling s the cantilever surface, 
L m its length and k its thermal conductivity, from our 
approximations and Fourier's law we obtain a linear tem- 
perature profile, 



T{y) = T 



an c (0)hu c y 



(39) 



The effective temperature to be used in Eq. (|27l) will thus 
be an average between T(0) and T(L m ), taking into ac- 
count that the mechanical mode profile is not distributed 
uniformly over y 



T 



T 



an c (0)huj c L r 



(40) 



where the adjustable parameter e is of the order of unity 
(e = 2 for an arithmetic mean). We include this effect in 
our estimates of x i n t ne next section. 



IV. EXISTING EXPERIMENTS ON 
PHOTOTHERMAL COOLING: ESTIMATING x 

In order to test our model and to gain some insight in 
the value of the key parameter x, we analyze the results 
from Ref. that, to the best of our knowledge, is the 
only publication reporting the observation of photother- 
mal cooling with enough experimental data to allow a 
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comparison with our theory. In Ref. [2^ the temperature 
of the cantilever mode is measured for different laser pow- 
ers, and for temperatures between 300 K and 32 K. All 
the material parameters but — can be retrieved directly 
from the aforementioned reference. 

Fitting the temperature for each different laser power 
allows us to fix an estimate for the value of x around the 
value 2 x 10~ 5 s m _1 . With these parameters, we obtain 
a phonon noise population of 1.4 x 10 4 , that is a phonon 
temperature of 5 mK. 

V. CONCLUSIONS 

In summary, we have extended the classical model of 
photothermal cooling by Metzger et al.Q to take into 
account the effect of quantum noise on cooling efficiency. 
Rather than build up a quantum model from a micro- 
scopic Hamiltonian, we treated the complex many-body 
thermal cantilever deformation effect as an effective time- 
retarded force. As one cools the mechanical motion to the 
ground state, noise due to both photon absorbtion and 
radiation pressure begins to play a role, and we explic- 
itly included these, and heating due to absorbtion, in our 



model. 

In the future it will be interesting to test the predictive 
power of our model against a larger experimental data-set 
than the one currently available. 
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